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Non-perturbative renormalization constants and light quark masses 
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We present the results of an extensive non-perturbative calculation of the renormalization constants of bilin- 
ear quark operators for the non-perturbatively 0(a)-improved Wilson action. The results are obtained at four 
values of the lattice coupling, by using the RI/MOM and the Ward identities methods. A new non-perturbative 
renormalization technique, which is based on the study of the lattice correlation functions at short distance in 
a;-space, is also numerically investigated. We then use our non-perturbative determination of the quark mass 
renormalization constants to compute the values of the strange and the average up/down quark masses. Af- 
ter performing an extrapolation to the continuum limit, we obtain mfS(2 GeV) = (106 ± 2 ± 8) MeV and 
mf^{2 GeV) = (4.4 ± 0.1 ± 0.4) MeV. 

Although the calculation of renormalization 
constants (RCs) can be performed in per- 
turbation theory, the implementation of non- 
perturbative renormalization (NPR) techniques 
has become an essential ingredient of lattice QCD 
calculations. These techniques allow to control 
the systematic error associated with the determi- 
nation of the RCs at the level of accuracy reached 
at present by the statistical and other systematic 
uncertainties in lattice calculations. 

In this talk we present the results of an ex- 
tensive NP calculation of the RCs Zv, Za, 
Zs, Zp and Zt of the bilinear quark operators, 
for the non-perturbatively C'(a)-improved Wilson 
action |^ . These results are obtained at four val- 
ues of the lattice coupling, in the range 6.0 < 
f3 < 6.45, by using the RI/MOM method [| and, 
for the scale independent RCs, the Ward identity 
(WI) approach g. A NP RI/MOM determina- 
tion of the RCs of the AF = 2 four-fermion oper- 
ators has been also presented by J. Reyes at this 
Conference [||. 



* Presented by Vittorio Lubicz at "Lattice 2002", Boston. 
tThis work is supported in part by the European Network 
"Hadron Phenomenology and Lattice QCD" , HPRN-CT- 
2000-00145. 



The quark mass renormalization constants 
computed in this study have been used to calcu- 
late the strange and the average up/down quark 
masses That calculation, which is completely 
0(a)-improved, uses NPR and involves an ex- 
trapolation to the continuum limit, represents at 
present one of the most accurate determination of 
the light quark masses within the quenched ap- 
proximation. 

Finally, we present the results of a preliminary 
numerical investigation of a NPR technique based 
on the study of the lattice correlation functions 
at short distance in x-space (XS) @]. 

1. RENORMALIZATION CONSTANTS 

The lattice parameters used in this study are 
summarized in table ^ The values of the inverse 
lattice spacing given in the table have been de- 
termined from the study of the static quark anti- 
quark potential fixing in input a~^{(3 ~ 6) ~ 
2.0(1) GeV. 

In order to study finite size effects, two inde- 
pendent simulations on different lattice sizes have 
been performed at P — 6.0. The smallest size 
corresponds approximately to the same physical 
volume used at the larger values of (3. From this 



2 



Table 1 

Summary of the lattice parameters used in this work. 



p 


a-^(GeV) 




# K« 


# Confs 


6.0 


2.00(10) 


W X 52 


4 


500 


6.0 


2.00(10) 


243 X 64 


3 


340 


6.2 


2.75(14) 


243 X 64 


4 


200 


6.4 


3.63(18) 


323 X 70 


4 


150 


6.45 


3.87(19) 


323 X 70 


4 


100 




p=6.2 




p=( 


5.2 
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Figure 1. The scale independent RCs Zv, Za and 
the ratio Zp/Zs (left) and the RGI combinations 
Z^°', Z§^^' and Z^^' (right) as a function of the 
renormalization scale, at /3 = 6.2. The solid lines 
represent the results of a constant fit to the data. 

analysis, we find that finite volume effects in the 
calculation are limited in the range between 1% 
(Zv, Za) and 3% (Zs, Zp, Zt). 

The RI/MOM determination of the RCs is 
based on the study of the lattice correlation 
functions at large and in the chiral limit. 
The renormalization condition, in the RI/MOM 
scheme, ensures that leading 0(a)-discretization 
effects vanish in these limits [|l|. Therefore, the 
improvement of the external quark fields and of 
the composite operators is not necessary for this 
determination. The on-shell improvement of the 
vector and axial-vector current operators is re- 
quired, instead, when the RCs are computed from 
the study of the lattice Wis. In this case, we 
use the coefficients cy and ca determined non- 
perturbatively in refs. (|J^. 

The results for the scale independent RCs Zy 
and Za, and for the ratio Zp/Zs, obtained with 
the RI/MOM method from the simulation at 
(3 = 6.2 are shown, as an example, in fig.|l| 
(left) as a function of the renormalization scale. 
The good quality of the plateau indicates that 
0(a^) discretization effects are well under con- 



• ZP subtracted 

• ZP'*'^' unsubtracted 
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Figure 2. Zp'^' as obtained with and without the 
subtraction of the Goldstone pole, at /3 = 6.4. 

trol. In figji (right) we also show the renor- 
malization group invariant (RGI) combinations 
Z§(^^ = Zo(^i)/Co{t^), for O = S,P,T, where 
the evolution function CoifJ.) is introduced to 
cancel the scale dependence of the RCs. In the 
RI/MOM scheme, these functions are known at 
the NLO in perturbation theory for Zt and at 
the N^LO for Zs and Zp From the quality 
of the plateau, discretization effects appear to be 
very small, even at quite large values of (a/z)^. 
Our estimates of the RCs have been obtained by 
fitting to a constant the results shown in fig.|l| 
(and similarly for the other (3s). The uncertainty 
due to the residual scale dependence observed in 
the plots has been used to evaluate the systematic 
errors. 

The RI/MOM approach relies on the fact that 
NP contributions to the Green functions vanish 
asymptotically at large p^. Special care, how- 
ever, must be taken in the study of the pseu- 
doscalar Green function T p since in this case, due 
to the coupling to the Goldstone boson, the lead- 
ing power suppressed contribution is divergent in 
the chiral limit In order to improve the 

convergence of Tp to its asymptotic behaviour 
at large p^, we have subtracted the Goldstone 
pole contribution by constructing the combina- 
tions IT2I 



SUB/ _ ^1 
P (nil, 1712) 



mirp(mi) - m2Tp(m2) 



(1) 



mi — TO2 

where mi and m2 are different quark masses. The 
effect of this subtraction on Zp'^^ is shown in 
fig.|[ Though both the unsubtracted and the sub- 
tracted determinations of Zp'^^ converge to the 
same value at large p^, a clear plateau is never 
reached in the first case, in the region of momenta 
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Table 2 

Values of the RCs obtained with the RI/MOM method. The scheme and scale dependent constants, Zs, Zp and 
Zt are given in the RI/MOM scheme at the scale fi = 1/a. 



a 


6.0 6.2 6.4 6.45 




6.0 6.2 6.4 6.45 


Zv 


0.766(2) 0.775(3) 0.795(2) 0.797(3) 


Zs 


0.668(13) 0.677(9) 0.696(6) 0.707(8) 


Za 


0.804(2) 0.809(3) 0.824(2) 0.825(4) 


Zp 


0.535(12) 0.564(4) 0.600(5) 0.612(9) 


Zp/Zs 


0.804(6) 0.831(3) 0.862(3) 0.867(3) 




0.833(2) 0.847(2) 0.867(6) 0.867(5) 



Table 3 

VaJues of the RCs obtained with the WI method. 





6.0 6.2 6.4 


Zv 
Za 
Zp/Zs 


0.774(4) 0.789(2) 0.804(2) 
0.841(17) 0.811(5) 0.843(10) 
0.845(20) 0.863(5) 0.911(10) 



explored in this study. 

Discretization effects can also be investigated 
by comparing the results for the RCs obtained at 
different values of the lattice spacing. Though the 
RCs depend on the lattice spacing (the UV cutoff 
in the lattice regularization) , their dependence on 
the renormalization scale is the same for all a, 
and it is only fixed by the anomalous dimension 
of the relevant operators. Therefore, the scale 
dependence should cancel in the ratio 

Z(a,/i)/Z(a',^) = i?(a,a'), (2) 

up to discretization effects. In fig.^we show as an 
example the results for the RCs Zy and Zp ob- 
tained at four values of /9 as a function of the 
renormalization scale. Each RC, Z{a',fj,), has 
been rescaled by the factor R{a,a') of eq. (||), 
where we have chosen as a reference scale the 
value of the lattice spacing at /3 = 6.4. In fig.|| 
we observe that the results obtained at the differ- 
ent values of /3 all lie on the same universal curve, 
thus confirming that discretization effects are rea- 
sonably small, within the statistical errors. The 
figure also shows that the renormalization scale 
dependence is in very good agreement with the 
one predicted by the anomalous dimensions of the 
relevant operators. 

Our results for the RCs of the bilinear quark 
operators determined with the RI/MOM and the 
WI methods are collected in tables || and respec- 
tively. These results are shown in fig.0, where 
they are also compared with the predictions of 
1-loop boosted perturbation theory (BPT). We 
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Figure 3. Zv and Zp at the four values of the lattice 
coupling as a function of the renormalization scale. 
The results have been rescaled by the factor R{a, a') 
defined in eq. The solid lines represent the scale 
dependence predicted by the corresponding anoma- 
lous dimensions. 










Zp(l/a) 


i 


0.92 0.94 0.9S 0.90 1.00 1.02 
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Z.j.(l/a) 





0.90 0.93 0.94 O.962 

90 



90 



Figure 4. VaJues of the RCs obtained from the 
RI/MOM method (Riled circles), the WI method 
(empty circles) and 1-loop BPT (solid lines). 

observe a good agreement between the RI/MOM 
and the WI determinations for the cases in which 
the latter method can be used, namely Zy, Za 
and Zp/Zs- There is also a general good agree- 
ment between NP results and BPT. A notable 
exception is the RC of the pseudoscalar density 
Zp, whose NP determinations differ from the per- 
turbative prediction by approximately 10-15%. 
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2. LIGHT QUARK MASSES 

An important application of the NPR is the lat- 
tice determination of the quark masses. We have 
computed the strange and the average up/down 
quark masses by using the standard vector and 
axial-vector WI methods, for which the relevant 
RCs are XjZg and ZaIZp respectively. We have 
used the same sets of gauge configurations and 
quark propagators used for the calculation of the 
RCs. The physical values of the light quark 
masses have been fixed by using the experimen- 
tal pion and kaon masses. Particular attention, 
in the calculation, has been dedicated to control 
and reduce the systematic uncertainties, within 
the quenched approximation. The final results 
have been obtained by performing an extrapola- 
tion to the continuum limit. All details of this 
calculation are discussed in ref . . Here we only 
quote our final results, which read: 

mf^(2 GeV) = (106 ± 2 ± 8) MeV 

mp(2 GeV) = (4.4 ± 0.1 ± 0.4) MeV . (3) 

3. THE X-SPACE METHOD 

In the XS method [Q one imposes the following 
renormalization condition 

Zo(xo/a) (O(x)O(O)) 1,2^,2 = (0(a:o)0(0))eo„t(4) 

where xq is a short distance scale. The symbol 
(. . .)cont denotes the Green function renormalized 
in a continuum scheme (the simplest choice is to 
fix this function to its value in the free theory). 

With respect to the RI/MOM method, the XS 
approach presents an important theoretical ad- 
vantage: the correlation functions in eq. (^ are i) 
gauge invariant and ii) not affected by the pres- 
ence of contact terms for xq ^ 0. This allows 
to neglect, in the renormalization procedure, the 
mixing with non gauge invariant operators and/or 
operators vanishing by the equation of motion. 

As a preliminary study of the XS method 
we show in fig.^ the results for the ra- 
tios (y(a;)y(0))/(A(x)A(0)) and (P(a:)P(0))/ 
(S'(a;)5'(0)), as a function of the distance 
(points corresponding to the same value of 
have been averaged together). These quantities 




(x/a) 



(x/a)' 



Figure 5. Ratios of 2-point correlation functions as a 
function of the distance (x/a)^. For comparison, the 
RI/MOM estimates of {ZA/Zvf and {Zs/Zpf are 
also shown with solid lines. 

provide directly the scheme and scale indepen- 
dent ratios (Za/Zv)"^ and (Zs/Zp)"^. For com- 
parison, the RI/MOM estimates of this RCs are 
also shown in the figure. We note that, though 
the results obtained with the XS method are con- 
sistent with the expectations, they are affected 
by large systematic uncertainties. Understand- 
ing these uncertainties requires further investiga- 
tions. 
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